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Abstract. - We theoretically demonstrate how the resonant coupling between a pair of cavity- 
polaritons and a biexciton state can lead to a large single-photon Kerr nonlinearity in a semi- 
conductor solid-state system. A fully analytical model of the scattering process between a pair 
of cavity-polaritons is developed, which explicitly includes the biexcitonic intermediate state. A 
dramatic enhancement of the polariton-polariton interactions is predicted in the vicinity of the 
biexciton Feshbach resonance. Application to the generation of non-classical light from polariton 
dots is discussed. 



A number of quantum optical applications crucially rely on having a strong value of the 
effective photon-photon interaction mediated by the optical nonlinearity of the atomic or 
solid-state medium. The most interesting quantum physics appears in fact when the presence 
of a single photon is sufficient to significantly modulate the response of a device. As a 
simplest example, a stream of strongly antibunched photons is emitted by a cavity as soon as 
the nonlinear shift of the resonance by a single photon is larger than the cavity lincwidth |1I3) . 
This physics is even more intriguing when more complex devices are considered: photons 
have been predicted to fermionize into Tonks- Girardeau gases as soon as the impenetrability 
condition is satisfied in a one-dimensional geometry. To this end, both strongly nonlinear 
optical fibers [3] and arrays of many coupled nonlinear cavities [5] have been considered. 
Quantum phase transitions between a coherent phase of "superfluid" photons to a Mott 
insulator one have also attracted a significant interest [6]. So far, the main obstacle against 
an experimental realization of all these phenomena has been the lack of scalable optical 
devices with sufficiently large nonlinearities and weak enough losses. 

A great deal of the recent advances in the field of strongly correlated atomic gases have 
been made possible by the discovery of the so-called Feshbach resonance effect in atom- 
atom collisions [7]: the scattering cross section is dramatically enhanced when the energy 
of the pair of colliding atoms is resonant with a long-lived quasi-bound molecular state. In 
typical experiments, an external magnetic field is used to tune the energy of the quasi-bound 
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molecular state close to the energy zero of scattering states. At this point, the low-energy 
scattering amplitude diverges and a dilute ultra-cold atomic gas has the strongest possible 
interaction. Remarkable recent experiments with atomic Fermi gases have exploited such 
Feshbach resonances to demonstrate superfluid behavior at temperatures comparable to the 
Fermi temperature [3]. 

In this Letter, we show how one could use a solid-state analog of the Feshbach resonances 
to dramatically enhance the strength of the interactions between cavity-polaritons in a 
planar microcavity in the strong coupling regime [10) . Inserting realistic parameters into our 
analytic model, we estimate that the enhanced nonlinearity should enable the realization of a 
polariton dot - an anharmonic quantum emitter that does not require quantum confinement 
of carriers. Since they are not subject to the fluctuations in the confinement length that are 
responsible for the typical inhomogeneous broadening of quantum dot samples, polariton 
dots appear as ideal candidates to realize sizeable arrays of identical nonlinear cavities. 

During the last two decades, a significant literature has appeared on the physics of biex- 
citons in semiconductor materials, i.e. two electron-two hole bound complexes analogous to 
hydrogen molecules [5]. A recent work by M. Wouters [IT] has pointed out the possibility of 
exploiting biexciton Feshbach resonance effects to enhance the optical nonlinearity of planar 
microcavities: by adjusting the cavity length it is in fact possible to bring the energy of a 
two-polariton state in resonance with a biexciton state. Our work follows up on this pio- 
neering proposal: we develop a simple analytical model of the biexciton Feshbach resonance 
and we analyze the nonlinear and quantum optical phenomenology of specific semiconduc- 
tor devices. The theoretical model is inspired by the atom-molecule approach to strongly 
interacting degenerate quantum gases |12 [ ll3 j and provides analytical formulas for both 
the energy-dependent scattering amplitude in a planar geometry and the single-polariton 
nonlinearity in confined geometries. 

The model. — We describe the dynamics of the two-dimensional interacting bosonic 
fields describing the cavity-photon fc,ir(r), exciton ^>x,a{ Y ) and biexciton ^ s(r) in a planar 
cavity with an embedded quantum well by a model Hamiltonian H = Hxc + Hbxc in 
second quantized form. As in the atom- molecule case (T2HB], separating out the bound 
(biexciton) state is legitimate as long as its spatial extension (~ lOnm) is much smaller than 
the characteristic length over which the dynamics of the exciton and cavity-photon fields 
is taking place (~ l^m). The single particle dynamics of cavity photons and excitons is 
described by the term: 




where E£ x B arc the rest energies of a cavity-photon, an exciton and a biexciton, respec- 
tively. Mc,x,b are the corresponding effective masses within a parabolic dispersion approx- 
imation. In the typical experimental systems, the exciton and biexciton masses Mx,b are 
comparable to the free electron mass (mo) and are much larger than that of the cavity- 
photon (Mc ~ 10~ 5 m ). The spin variable a =t, I identifies the polarization states of the 
cavity photon and of the exciton. The biexciton state in contrast is a singlet. The dark 
exciton states are not included in the model since they do not contribute to the physics we 
consider. Vm,x),aa' ( r ) denotes the - possibly spin-dependent [25] - confinement potential 
acting on the cavity-photon and the exciton, respectively. SIr is the Rabi frequency corre- 
sponding to the exciton-photon dipole coupling. A typical example of the momentum-space 
dispersion of the coupled exciton and cavity-photons in a spatially homogeneous geometry 
Vx,c — is shown in FigfJJa): the exciton and cavity photon (dotted lines) are mixed into 
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the lower- (LP) and upper-polariton (UP) states (solid lines). On the scale of the figure, 
the dispersion of excitons is almost flat (dashed line). 

The principal process leading to the formation of a biexciton involves the absorption of 
a cavity-photon by an existing exciton and can be modelled by the following Hamiltonian 
®: 

Hbxc = ~ I d2Rd2r [g(r) & B pL)9xAR ~ r/2)* c ,_ ff (i? + r/2) + h.c] . (2) 

Even without a detailed knowledge of the microscopic structure of the biexciton, general 
physical arguments can be invoked to obtain important information on the coupling ampli- 
tude g(r): the optical process leading to the formation of a biexciton is limited to a spatial 
region around the existing exciton with a characteristic size on the order of the spatial size 
of the biexciton. Its value scales as the amplitude of the wave-function describing the 
relative motion of two excitons forming a biexciton a~j^~h{r/a,B)\ here, h(s) is determined by 
the attractive exciton-exciton interaction potential and has a peak value and a spatial exten- 
sion of order unity. The local value of the exciton creation matrix element is £l B . These two 
general observations lead us to the following generic form for the exciton-photon-biexciton 
coupling: 

9(r) = ^h(r/a B ). (3) 
a B 

If the exciton-photon Rabi frequency fl p is comparable to or larger than the biexciton 
binding energy, we can tune the cavity in a way to ensure that the upper-polariton branch 
is far enough away in energy to be safely neglected and, at the same time, that the biexciton 
resonance is comparable to twice the energy of polariton states near the bottom of the 
lower-polariton dispersion [Figflja,c)]. We can then rewrite the exciton-photon-biexciton 
coupling term H B xc m terms of the lower polariton field only: 



Hbxc = ~ V f d 2 R d 2 r ^ u£ P u c LP x 
2 _ 1 J a B 

x / ! ( r /a B )^(R)* LPi „(R-r/2)^_ J (R + r/2)+h.c. (4) 



X c 

Here, u LP are the Hopfield coefficients quantifying the exciton and cavity-photon com- 
ponent of the lower polariton; as we are focussing our interest to the lowest part of the 
lower-polariton dispersion only, the dependence of u^' p (k) on momentum k is here ne- 
glected and their value at k — is considered. The form @ of the polariton-biexciton 
coupling emphasizes the crucial fact that the final product of the biexciton disintegration 
consists of a pair of lower polaritons rather than an exciton and a photon |14j . 

In general, there are additional exciton scattering processes that are independent of the 
biexciton resonance, analogous to the so-called background scattering length in the atomic 
case [7]. As is the case in the theory of cold quantum gases, the effective potential describing 
such scattering processes does not need to be the physical exciton-exciton potential. Any 
model potential can be chosen, provided it correctly reproduces the scattering physics in 
the region of parameters of interest once all direct and exchange processes are taken into 
account [15]. For the sake of simplicity, in the following we shall choose a repulsive, short- 
ranged form for the exciton-exciton potential t><x,CT'(r — r') = v a ,a> (5 2 (r — r'). If the relevant 
physics is limited to the lower-polariton branch, the same reasoning underlying Q leads to 
a local polariton-polariton potential of rescaled intensity v a ^> \u^ P \ 4 . 

Polariton-polariton scattering amplitude. Since our principal focus is on the 
biexciton-Feshbach resonance enhancement of the polariton-polariton interactions, we first 
neglect the contribution from the direct exciton-exciton scattering. Straightforward manipu- 
lations can be used to resum the diagrams contributing to the polariton-polariton scattering 
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amplitude stemming from Eq. (2). Approximating the polariton dispersion with a parabola 
of mass Mlp, the scattering T- matrix element in the singlet channel has the resonant form: 

TU{E) = E-E° B + M^ jip \og[E max /E] + zMjf \g\* (5) 

This expression for the T-matrix element was first discussed in |16j in the context of atom- 
atom scattering in two-dimensional geometries. The coefficient g is here defined as g — 
hilpapUppUpph with h — a~^ 2 J d 2 r h(r /as) is of order unity. No scattering occurs instead 
in the triplet channel. 

From equation ([5]) , it is easy to recognize several features typical of a Feshbach resonance 
in two-dimensional scattering. The radiative correction to the resonance energy has a log- 
arithmic dependence on the energy of the form A ra d = Mpp \g\ 2 /(AttH 2 ) \og[E max / E]; the 
high-energy contribution to A rac i depends very much on the details of the system and has 
been absorbed into the ultraviolet cut-off E max . The imaginary part of the denominator 
accounts for the radiative decay of the biexciton into a pair of propagating polaritons with 
an energy-independent linewidth T ra d = Mpp \g\ 2 /2h 2 comparable to the radiative shift 
A. rac i- Inserting realistic parameters from GaAs-based structures, namely Kip ~ 5 meV 
and ap — 20 nm, the resulting radiative linewidth T ra d into a pair of polaritons is on the 
order of a few /zeV, somehow smaller than the value quoted in [17] . Additional contributions 
to the biexciton decay rate could be included in an effective T to t — ^rad + ^add- A significant 
contribution to the additional decay rate T a dd may stem from dissipative coupling of the 
biexciton to bulk photons and/or interface polaritons. For simple planar geometries, we can 
estimate this contribution to be on the order of 10 fieV. 

At the mean-field level the dynamics of the coupled cavity-photon, exciton and biexciton 
quantum fields can be described in terms of classical fields &c,x,b — (^c,x,b(^)} ■ This 
approach has been already proven to be useful in studies of biexciton optics [HUH] as well 
as to describe the dynamics of atom-molecule conversion in ultracold atomic gases [12) 115]. 
In contrast, we shall now focus on a spatially confined system where quantum fluctuations 
are instead strong and the mean-field approximation fails. Our goal is to demonstrate that 
the strong polariton-polariton interactions originating from the biexciton Feshbach reso- 
nance can translate into a significant polariton blockade effect [3] with peculiar antibunching 
features in the transmitted light. 



Zero-dimensional polariton dots. — Zero-dimensional polariton dots [19 21 arc 
characterized by a strong spatial confinement of the polariton wavefunction in all three di- 
mensions, which gives rise to a series of discrete polariton states. In typical experimental 
set-ups, the in-plane confinement of polaritons is obtained by suitably designing a photonic 
structure whose discrete photonic eigenmodes have a small transverse mode profile. Con- 
finement of the bare exciton or the biexciton is instead not necessary as it gets inherited 
from the photonic one. 

If the confined polariton modes are well separated in energy as compared to their 
linewidth and all nonlinear couplings, we can restrict our attention to a single polariton 
mode. This assumption requires that we limit ourselves to geometries where the photon 
is confined on a length scale of the order of a /im. A moderate spatial anisotropy of the 
confinement is enough to have a sufficient spectral separation (> lmeV) of the two lowest 
linearly-polarized polariton modes [25] , so that the spinor part of the (spinorial) wave- 
function (f>pp(r,cr) becomes trivial 0lp ( r ,t) = <f>Lp{r,l) — ^75 </>z,p(r). 

The creation operator of the single confined linearly polarized mode will be denoted 
by p. Nonlinear processes occur via resonant scattering onto the intermediate biexciton 
state. Given the form ^ of the polariton-biexciton coupling and the very large value of 
the biexciton mass as compared to the polariton one, a single biexciton state of spatial 
wavefunction 4>b(j) proportional to the square 4>pp(v) 2 will be involved in the dynamics. 
The corresponding biexciton creation operator will be denoted as b. We emphasize that the 
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biexciton states are not confined in our model; the wave- function <f>B (r) is determined by the 
collective coupling amplitude to the relevant biexciton eigenmodes contributing to confined 
polariton-polariton scattering. 

Limiting our attention to these polariton and biexciton states, the resulting Hamiltonian 
can be written in the concise form: 

H = E LP p^p + E B ~b r b+^Vp^p^pp + G tfpp + p^p r b . (6) 

Here, Elp and Eb are the energies of the confined lower-polariton and biexciton states, 
respectively. The nonlinear interaction coefficients V and G are evaluated by inserting the 
lower-polariton <f>Lp(r,o~) and biexciton </>p (r) wave- functions into the interaction Hamilto- 
nian and neglecting contributions from all other states, 

V = \Y, v ^'\ u lp\ 4 J^\Mr)\ 4 (7) 

cr.tx' 

B -uf P u c LP J d 2 Rd 2 r/ l (r/a B )^(R)0 iP (R-r/2,t)0Lp(R + r/2,|) 

(8) 



G = 



as 



-KlRdBufpU^p h 



1/2 

d 2 r|^p(r)| 4 



An important direct consequence of ([7]) and (|8|) is that the background interaction energy 
decreases as l/£ 2 for increasing confinement size £, while the polariton-biexciton coupling 
has a slower scaling proportional to l/£ only. 

Previous calculations of the background exciton-exciton scattering [15] suggest that the 
background interaction energy v a is of the order of R*a 2 x , where R* is the exciton Rydberg 
energy (on the order of 10 meV for standard GaAs quantum wells) and ax is the exciton 
Bohr radius (on the order of 10 nm). As the radius ap ~ 20nm of Wannier biexcitons is of 
the same order as the exciton radius ax and the Rabi energy Sip sw 2meV is typically of the 
same order as the exciton Rydberg R* , the resonant contribution G is expected to dominate 
over the background scattering contribution V as soon as the polariton dot is much larger 
than the physical size of excitons and biexcitons I ^> ax,B- As the confinement length I of 
typical polariton dots is on the order of 1 /xm, this latter condition is generally very well 
satisfied. 

Inserting these values of the parameters in ([8]), one obtains a value of G on the order of 
30 /ieV, which is extremely promising in view of applications: confined polariton states with 
a linewidth of 70/ieV have been recently demonstrated [3T]. On the other hand, this value 
of the nonlinear coupling is safely below the energy distance from the neighboring polariton 
states in the dot; the typical energy splitting between confined polariton states in a box 
potential with a transverse size of the order of 1/mi are in fact on the order of a few meV, 
which guarantees a posteriori the validity of our single-mode calculations. As compared to 
the original proposal [3] , the use of the biexciton Feshbach resonance allows us to relax the 
constraint on the size of the polariton dot. 

Few-body states. — As the Hamiltonian ([6]) conserves the total number of excitation 
Not = Nlp + 2Np, it can be easily diagonalized within each N to t manifold, giving a 
ladder of mixed levels. An example of this level ladder is shown in FigQJb) for the most 
remarkable case, 2E^p + V = Eb, where the biexciton is on resonance with the (slightly 
shifted) polariton states. The Ntot = manifold consists of the vacuum state |vac) only, 
while the Ntot = 1 manifold consists of a single one-polariton state \p) with energy £pp. 
The Ntot = 2 manifold contains of a doublet of eigenstates which are the symmetric and 
antisymmetric combinations of a biexciton \B) and a pair of polaritons with identical linear 
polarization, \pp). The energies of these states are symmetrically split by an amount ±\2G, 
typically much larger than V. In this way, the harmonicity of the level ladder is broken: the 
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Fig. 1: Upper-Left (a) panel. Solid lines: dispersion of the polariton frequency as a function 
of in-plane momentum; the exciton-cavity detuning is chosen such that the lower polariton is on 
(Feshbach) resonance with the biexciton energy (horizontal dotted line). Dashed lines: bare cavity 
photon and exciton dispersion. Upper-Right (b) panel: schematic diagram of the states belonging 
to the Ntot = 0, 1,2 manifolds exactly on Feshbach resonance 2Elp + V = Eb- Lower-Left (c) 
panel, k = anticrossing of the lower and upper polariton branches as a function of the bare 
cavity frequency Eq. The horizontal dashed line indicates half the biexciton energy; the Feshbach 
resonance is located at its crossing point with the lower polariton branch. Lower-Right (d) panel. 
Position of the different states belonging to the N to t = 2 manifold as a function of the bare cavity 
frequency Eq. Dashed line: bare \pp) and \B) states. Solid lines: states originating from the 
mixing of the biexciton state \B) and the \pp) state. The vertical dotted line indicates the Feshbach 
resonance point. Parameters: Ex = 1-4 meV; £Ir = 1.5 meV; G = 0.03 meV; V — 0. 



transitions from the Ntot — 2 to the Ntot = 1 manifold occur in fact at energies -Elp ± \/2G 
different from the energy E^p of the Ntot = 1 to N tot = transition. 

These simple considerations straightforwardly extend to the eigenstates forming the 
higher Ntot > 2 manifolds, which then consist of linear superpositions of states with N to t — 2n 
polaritons and n biexcitons (0 < n < Not/ '2) with suitable weights. In particular, it is in- 
teresting to note that, differently from the atomic case, the absence of deep molecular states 
other than the biexciton state should suppress all those three- and many-body recombination 
processes that generally limit the lifetime of strongly interacting atomic Bose gases [2"1] . 
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Far from resonance \2Elp — Eb \ G, the polariton and biexciton states are no longer 
efficiently mixed by the polariton-biexciton coupling and their energies approach their un- 
coupled values. In this regime, the harmonicity of the level ladder is only broken by a much 
weaker background scattering term of magnitude V <C G. A complete picture of the an- 
ticrossing between the biexciton and polariton states in the N to t = 2 manifold is shown in 
FigQJd): while in the atomic case the difference in magnetic moment of the molecular and 
atomic states allows to tune their relative energy by means of an external magnetic field, 
here the resonance between the polaritons and the biexciton can be obtained by suitably 
adjusting the cavity length. As it is shown in FigQJc,d), the polariton energy strongly de- 
pends on the cavity-photon energy (and in turn on the cavity length), while the biexciton 
energy only depends on the microscopic structure of the quantum well. 

Optical signatures of polariton blockade. — The pronounced anharmonicity of 
level scheme of FigQJb) has direct implications for the optical properties of the polariton 
dot. For the sake of simplicity we restrict here our attention to the G 3> I\ T a dd limit, 
where T is the radiative linewidth of the polariton states and T a dd is the additional, non- 
radiative contribution to the biexciton linewidth. The driving laser is assumed to be fully 
coherent, weak and resonant with the transition to the one polariton state \p), which then 
gets effectively populated. On the other hand, the Feshbach-resonance-induced nonlinearity 
shifts both states in the N to t = 2 manifold far from resonance with the laser by an amount 
\plG larger than the linewidth of the relevant states. As a consequence of this polariton 
blockade effect, the Ntot — 2 manifold remains almost unpopulated: the first polariton is 
able to effectively block the entrance of the second. In a transmission geometry where only 
light emitted by the polariton dot is collected, this polariton blockade effect is clearly visible 
as a complete anti-correlation between transmitted photon^. 

The decay rate Y of the lowest energy cavity-polaritons is determined exclusively by 
the cavity-losses. While the lowest value observed to date is V ~ 70/ieV, there are no 
physical constraints that prohibit reaching much smaller values by increasing the cavity 
quality factor. On the other hand, the intrinsic biexciton decay rate Y a dd into bulk photons 
or surface polaritons can be estimated to be on the order of 10/ieV. Unless some specific 
effort is done to suppress it in a suitably designed photonic structure, this biexciton decay 
rate is most likely to set the ultimate limit to the antibunching that is achievable in polariton 
dots. 

A recent experiment [22] has demonstrated how dissipative collisional processes offer 
a viable alternative route to create strongly correlated atomic gases. The possibility of 
extending this idea to the photonic case was discussed in [33] . In our system, the use of such 
dissipative nonlinearities provides a simple way to obtain a strongly antibunched emission 
in the opposite regime T a dd G ^> T as well: while the energy spacing between the levels 
remains almost constant, the harmonicity of the level ladder in this case is broken by the 
strongly varying dissipation rate. For instance, the strong biexciton damping T a dd allows 
the decay rate of the two-polariton \pp) state T2 — 2T + 2G 2 /T a dd to be much larger than 
the decay rate T of the one-polariton \p) state. As a result, the excitation of the system 
under a weak and coherent pump, is limited to the Ntot = 1 manifold, which produces a 
stream of antibunched photons in the emission. 

Conclusions. — We have theoretically investigated the physics of coupled exciton, 
photon and biexciton quantum fields in planar semiconductor microcavities. The Feshbach 
resonance on the biexciton state can be exploited to dramatically enhance the amplitude of 
polariton scattering processes. A strong photon antibunching is predicted for the emitted 
light from strongly confined polariton dots. If our predictions are verified experimentally, 

J If the polarization degeneracy is not lifted, the Feshbach enhancement of the polariton-polariton in- 
teractions is restricted to the singlet channel, which limits the blockade effect to polaritons of opposite 
polarizations. This dramatically suppresses the anti-bunching features, but leaves interesting polarization 
correlations in the emitted light. 
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biexciton physics will become the cornerstone for a new generation of opto-electronic devices 
working at the single-photon level, with a number of applications from fundamental science 
to quantum information processing. 

* * * 

The core of this work was performed while I.C. was visiting the QPG group at ETH- 
Ziirich. Continuous exchanges with Michiel Wouters, Cristiano Ciuti, Giuseppe La Rocca, 
Andreas Reinhard, Hakan Tiireci, Antonio Badolato, Dario Gerace, and Rosario Fazio are 
warmly acknowledged. 
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